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Abstract 

Is studied asymptotic expansion for solution of singularly perturbed equation for 
Markov random evolution in IR d . The views of regular and singular parts of solution 
are found. 
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1 Introduction 

A Markov random evolution (MRE) is created by a solution of the evolutionary 
equation in Euclidean space M. d , d > 1 

du £ (t)/dt = v(u E (t);3s(t/e)) 

with the ergodic Markov switching process se(t),t > on the standard (Polish) 
phase-space (E,£) by the operator Q(x,B),x G E, B G £ that defines transition 
probabilities of a Markov chain se n , n > 

Q(x, B) = P{as n+ i G B\si n = x}. 

The operator of transition probabilities Q is defined by 

Qf(x) = J Q(x,dy)f(y),xeE, (1) 
for any bounded measurable real valued / defined on E. 
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We will see later that the equation for the regular and the singular parts of a 
random evolution are defined by the generator (1) of a uniformly argodic Markov 
switching process. The Banach space 13(E) is splitted onto the two subspaces [7]: 

B(E) = N Q @R Q , 

where Nq := {ip : Qip = 0} is the null-space of Q, and Rq := {ip : Qp = Tp} is the 
range of Q. 

We define the projector II : N Q := UB{E),R Q := (I -U)B(E); U(p(x) := (pl,ip:= 
j E p(x)ir(dx), where the stationary distribution it(B),B G £ of the Markov process 
dd(t),t > satisfies the relations [1] 

n(dx) = p(dx)m\(x) /m, 
m = I mi(x)p(dx). 
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p(B),B £ £ is the stationary distribution of the Markov chain ae n ,n > 0, given by 
the equation 

p{B) = j Q(x 1 B)p(dx),p(E) = l. 

E 

Let us consider the Banach space B(M. d ) of real- valued test-functions p(u),u £ M. d 
which are bounded with all their derivatives equipped with sup-norm 

1 1^| | := sup \<p(u)\ < C v . 



The random evolution in £>(R ) is given by the relation 

$f(u,z) := E[ip{u e (t))\u £ (0) = n,aa e (0) = x]. (2) 

The asymptotic behavior of MRE (2) as e — > is investigated under the assump- 
tion of uniformly ergodicity of the Markov switching process as(i) described above 
and under the assumption of the existence of a global solution of the deterministic 
equations 

du x (t)/dt = v(u x (t); x), x € E. 
Let us consider the deterministic evolution 

$ x (t,u) = ip(u x (t)),u x (0) = u. 

It generates a corresponding semigroup 

Y t (x)ip(u) := <p(u x (t)),u x (0) = n, 

and its generator has the form: 

d 

Y(x)ip(u) = v(u;x)p'{u) :=^2v k (u;x)(p' k (u), 

k=l 
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tp' k {u) := dtp(u)/du k ,<p(u) e C°°(M d ). 
By the average principle [TO] the weak convergence 

u £ (t) => u(t),e -> (3) 

takes place. The average limit evolution u(t),t > is defined by a solution of the 
average equation 

du{t)/dt = v(u(t)). 
The average velocity v{u) , u £ M. d is defined by 

v(u) = / v (u; x)Tr(dx) 



(i.e. by the average of the initial velocity v(u;x) over the stationary distribution 
n{B),Be £). 

The rate of convergence in (3) can be investigated in two directions: 
i) asymptotic analysis of the fluctuations 

C(t) = u %t)-u(ty, (4) 



ii) asymptotic analysis of the average deterministic evolution (2). 

The asymptotic analysis of fluctuations (4) leads to the diffusion approximation 
of the random evolution [TO] . 

The asymptotic analysis of evolution (2) is realized in what follows by constructing 
the asymptotic expansion in power of the small parameter series e — > 0(e > 0) in the 
following form (r = tje): 

oo 

® £ t (u,x) =/ ) (i) + ^£ k [/ ) (f)+/ ) (T)]. (5) 
k=l 

The asymptotic expansion (5) contains two parts: 

i) the regular term u £ (t) := «(°)(t) + *£T>1 £ k u {k) (t), 

ii) the singular term (boundary layer) w £ (t) := Ylk*>i ^ k w^ k \r), t = tje. 

In addition the initial condition: 

J°)(0) = <p(u)l 

has to be valid for any x G E,u G 

It's well-known (see, e.g. [8]), that the evolution, determined by a test-function 
<p{u) € C°°{R d ) (here <p{u) is integrable on W 1 ): satisfy the system of Kolmogorov 
backward differential equations: 

^t(u,x) = [e- 1 Q + Y]m(u,x), 
$ £ (u,x) =<p(u). 
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Asymptotic expansions with 'boundary layers' were studied by many authors 
(see [21 |3l[12])- In particular, functionals of Markov and semi-Markov processes are 
investigated from this point of view in [61 [9l lllj. 

In this work we study system (6) with the first order singularity. To find asymp- 
totic expansion of the solution of (2) we use the method proposed in [3l I12j . The 
solution consists of two parts - regular terms and singular terms - which are de- 
termined by different equations. Asymptotic expansion lets not only determine the 
terms of asymptotic, but to see the velocity of convergence in hydrodynamic limit. 

Besides, when studying this problem, we improved the algorithm of asymptotic 
expansion. Partially, the initial conditions for the regular terms of asymptotic are 
determined without the use of singular terms, i.e. the regular part of the solution may 
be found by a separate recursive algorithm; scalar part of the regular term is found 
and without the use of singular terms. These and other improves of the algorithm 
are pointed later. 

2 Asymptotic expansion of the solution 

Let P(t) = e Qt = {Pij{t);i,j G E}. Put ttj = lim Pij (t) and -R = {JT(Pij(t) - 

t — >oo 

7Tj)dt;i,j G E} = {ry;i, j G E}. 

Let II be a projecting operator on the null-space Nq of the operator Q. For any 
vector g we have Tig = gl, where g = (g, n), 1 = (1, . . . , 1). Then for the operator Q 
the following correlations are true (see [7j, chapter 3) 

UQU = 0, 

QR = R Q = IL-I. 

We put: 

ex Po (Qt) := e Qt -U, 

Theorem 2.1. The solution of equation (6) with initial condition <&q(u,x) = y>{u), 
where <p{u) G C°°{R d ) and integrable on W 1 has asymptotic expansion 

oo 

x) = u {0) (t) + Y, £n ( u(n) (*) + w(n) • (7) 

71=1 

Regular terms of the expansion are: u^(x,t) - the solution of equation 

^u {0 \t)-UYUu^{t) = (8) 
with initial condition u^°'(0) = <p(u), 

+ c (1) (t) := R JLuW{t)+c (1) {t), 
fork>2: 

u ( k )(t) = J R Ln (fc - 1) (t)+c (A;) (t) 



u 



(i) 



(t) = i? 



s» (0, w 



W°>(i) 
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where c^ k \t) G Nq, 

J» M -#)(v-'( t+ v(.)),o )+ / g^glM *-/ ^* > o, 

here V(u) = J V" —1 (u>) is i/te backward function for V (u) , 

k — 1 ' c— ' jk—i—n 

AM) = E E^- 1 )^^ " i " n + ^nv^n^— - c «( t ). 

i=0 n=l 

77ie singular terms of the expansion have the view: 

w^'(t) = expo(QT)Ytp(u), 

for k > 1 : 

r 

„«> W ^(Q^>(0) + /^ W r-.))V^(.)d.- 



oo 

n J Yw { ~ k ~ 1 \s)ds. 

T 

Initial conditions: 

C C°)(0) = p(«), 
^(^(O) = -UoL^ftt), 
c (1) (0) = 0, 

/or > 1 : 

w {k) (0) = -RolLu^iO), 
c ( fc )(0) = Vt? (fc - 1} (0), 

iw/iere ^(O) = -i? V^(n), 

#'(0) = floLu^- 1 )^) + J R V^ (fc - 1) (0)+ 

nV(55( fe - 1 )(A))' A | A=0 , 
(^ (fe) (A))' A | A= o = flgLu^CO) + j r2q 1 ^-i)( )+ 

Remark 2.2. The initial conditions for the regular terms of asymptotic are determined 
without the use of singular terms, i.e. the regular part of the solution may be found 
by a separate recursive algorithm (comp. with [6]). 
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Proof of Theorem 2.1: Let us substitute the solution §f(u,x) in the form (7) to 
the equation (6) and equal the terms at e degrees. We'll have the system for the 
regular terms of asymptotic: 



'q«(°) = 

Q u (k) = f t u {k -^ - Yu^-V := k^ 1 ', k > 1 



(9) 



and for the singular terms 

dw £ dw £ dr _ 1 dw £ 



dt dr dt 



= e 



dr 



(e- 1 Q + W)w £ 



Thus, from ^f- = (Q + eY)w £ we obtain: 



£wW - QwW = Yw^-V, k > 1. 



(10) 



From (9) we have: it (0) (t) G Nq. 

The solvability condition for uW (t) has the view: 

uQUu (1 \t) = o = |-« (0) m - nvn«(°)(t). 

at 

So, we have equation (8) for u^°\t). 
For (t) we have: 

« (1) (t) =i? Ln (0) (t) + c (1) (t). 
Using the second equation from (9) we obtain: 

u {k \t) = i? L« (fe-1) (*) + c (fc) (t), 

where c^(t) G TVq. 

To find C ( fc )(t) we'll use the fact that u^\t) G Nq. Let us put c^(t) = «(°)(t). 
For the equation 

Q„(2)(t) = _ Vu (i)( t ) = li? [4c (0) (i) " Vc(°)(t)l + 

Ot li .li 



dt dt 



dt 



Vi? 

we use the solvability condition 



d 



Ji c (0) (t )_ Vc (°)(t) 



VcW(t) 



nQiW 2 )(t) = o = -^(t) - Yc {l \t) + ni? n^ c (°)(t) -n/2oVn^- c (°)(t)- 

Ci6 Ci/t- til/ 



We find: 



nVi? n^c(°)(t) + UYRoYUc {0) (t). 

j t c [1 \t) - Yc {1 \t) = -nvi? vn c (0) (t). 
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By induction: 

, k— 1 k—i jk—i—n 

j/ k \t) - Vc<*>(t) = £(-!)*(*: l)nVfl V"n^— -c«(t), fc > 0. 

j=0 n=l 

So, we have the following equation for c^ k \u,t): 



here 



j=0 n=l 

To find a solution we should write down a system 

(ft du dc( fc ) 



1 i;(u) C k (u,t)' 
The independent integrals of this system are: 



i> (it) 

cW(M)+ m^du=c 2 . 

J v[u) 

As soon as c( fc )(it, t) is only in one of the first integrals, we may present the solution 
in the form: 

c^u,t) = f k (t + f^)- f^du,k>0, 

V J v ( u )J J v ( u ) 

where f k is any differentiable function. Using initial condition for (u, t) we find a 
condition for f k : 

A (/*) =<*><„, 0)+ f£i(^l du . 

\J v{u)J J v(u) 

We may put now V(u) = f and make a change of variables w = V(u). So, 
it = V~ 1 (w) and we have: 

fk ( W ) = cM(v-\w),o)+ r^(y-^o)dw 



v(V 1 (w)) w 
Thus, we obtain 
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Initial conditions for c^(u,0) are found later through Laplace transform for the 
singular terms of asymptotic. 

For the singular terms we have from (10): 

W W( T ) = exp (QT)w {l) (0). 

Here we should note that the ordinary solution iuW(r) = exp(Qr)w^ (0) is cor- 
rected by the term — Eht;W(0) in order to receive the following lim w^Or) = 0. We 

choose this limit to be equal for all singular terms, that may done due to uniform 
ergodicity of switching Markovian process. 

The following statements are made using a method proposed in [2j . For the second 
equation of the system the corresponding solution should be 

r 

-«W=« W (^)(0) + /« W (g(r-.))V^-) W *, 



where the homogenous part has the following solution 

w ( k \ T ) = exp o {Qr)w^{0). 

But here we should again correct the solution, in order to receive the limit 
lim w(*)(t) = 0, by the term -II f°° Ww^^ (s)ds. 
And so the solution is: 

r 

u.( fc )(r) = ex Po (QT)w {k) (0) + / exp (Q(T - s))Vu; (fc " 1) (s)ds- 



oo 

n J Vw^ k ~ 1 \s)ds. 



We should finally find the initial conditions for the regular and singular terms. 
We put c (°)(i) = «(°)(t), so c(°)(0) = «(°)(0) = f(u). 

From the initial condition for the solution u £ (0) = u^(0) = <p(u), we have to 
determine u^ k \0) +w^(0) = 0, k > 1. Let us rewrite this equation for the null-space 
Nq of matrix Q: 

Hu (fe) (0) + lW fc) (0) = 0,k > 1, (11) 
and the space of values Rq: 

(i - n)uM(o) + (/ - n) w ( k \o) = 0, k > 1. (12) 

So, for k = 1 we obtain: 

^(0) = J R L^°)(0) + c (1) (0) = 
= (/ - n)i2oLy)(u) + n c (1) (o), 



w W(o) = (i-il>«(o). 

Thus, cW(O) = OjitfWCO) = -i? Ly?(u). 
By analogy, for k > 1: 

u ( fc )(0) = iJoLu^*- 1 )^) + c (fc) (0) = 

= (J - n) J R Ln( fc - 1 )(o) + n c (fe) (o), 



(fc) 



(o) = (/ - il>W(o) - n / W* -1 ^)* 



o 



Functions u;^ 1 ^(s),u^' 1 ^(0) are known from the previous steps of induction. 
So, we've found lW fc )(0) in (11) and (I - IL>( fe )(0) in (12). 

Now we may use the correlations (11), (12) to find the unknown initial conditions: 

oo 

c ( fe )(0) = f Yw ( - k - 1 \s)ds, 







In [6] an analogical correlation was found for c^ k \0). To find c^ k \0) explicitly 
and without the use of singular terms we'll find Laplace transform for the singular 
term. The following lemma is true. 

Lemma 2.3. Laplace transform for the singular term of asymptotic expansion 

oo 

uj(*)(A) = J e- Xs w^{s)ds 



has the view: 

u>W(A) = (A - n + (r + nr^'HVMu)], 

w( k )(\) = (a - n + (R + n)- 1 )- 1 L U ( fc - 1 )(o)+ 
(A - n + (Rq + ny^vw^ix) + -nv^-^A) - ^-^(o)], 

A 



where 



w 



(fc) 



^W(O) = -R Vp(u), 
(w W WYx\x=o = -R 2 o^(u), 
(0) = flolLu^-^O) + J R V^ (fc " 1) (0)+ 

nv(^( fc - 1 )(A))' A | A=0 , 



(^ (fc) (A)) A | A=0 = ^L^" 1 )^) + J R2Q 1 ^ fc - 1 )(0)+ 
i? V(^( fc - 1 )(A)) A | A=0 . 



9 



Proof. 

oo oo 

{yW(A) = J e- Xs w {1) (s)ds = J e- Xs [e Qs - U]ds ■ w {1) (0) = 



= (A - n + (R + n)- 1 )- 1 ^^)], 

where the correlation for the resolvent was found in [7]. 

55W(0) = -R Yip(u), 
(w^(X))' x \ x=0 = Jim R ^\- Ro [-Y^(u)] = -R*V<p(u). 

A— >0 A 

For the next terms we have: 

w( k \\) = (A - n + (R + n)- 1 )- 1 W £ - 1 )(o)+ 

(A - n + (Rq + lT)~ 1 )~ 1 Yw( k ~~ 1 ') (A) + ^nV[w (fc - 1) (A) - w (fc ~ 1} (0)], 

A 

here the last term was found using the following correlation: 

oo oo oo 9 

J e- Xs J Yw ik - 1) (9)d9ds = J J e- Xs Yw ik - 1) (9)dsd9 = 

s 

\ {er xe - l)Yw^ k - 1 \9)dd = iv^-^A) - w (k - l \{))\. 



o 



So, 



w {k) (0) = Rohu^iO) + J R W fc - 1) (0) +IIV(^ (fc - 1) (A)) , A | A= 



o- 



(^ (fc) (A)) A | A=0 = Rlhu^iO) +j r2q 1 ^-D(o)+ 
R Y(w^ (A)) A j A=0 - Hm { ^UY[w^ (A) - (0)] - 

A^O ^ A z 

inv(^ fc - 1 )(A))' A 

where the last limit tends to 0. 
Lemma is proved. 

So, the obvious view of the initial condition for the c^ k \t) is: 

c ( fc )(0) = Vw (fe - 1} (0). 

Theorem is proved. 
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3 Estimate of the remainder 

Let function ip(u) in the definition of the functional $f belongs to Banach space of 
twice continuously differentiable by u functions C 2 (M. d ). 
Let us write (6) in the view 

&(t) = $ £ (t) - $|(t) (13) 

where = u^°\t) + e(u {1 \t) + w (1) (t)) + e 2 (u (2) (t) + w (2) (i)), and the explicit 

view of the functions u^(t),w^(t),i = 0,2, j = 1,2 is given in Theorem 2.1. 

By theorem 3.2.1 from [7] in Banach space C 2 {R d x £7) for the generator of 
Markovian evolution L £ = e~ l Q + V, exists bounded inverse operator (L e ) _1 . 

Let us substitute the function (13) into equation (6): 

— <l £ - L £ & = — $g - L £ <$>% := e(9 e . (14) 
dt dt 2 2 v ' 

Here e0 £ = e[f t u^ - eY(u^ + w^)]. 

The initial condition has the order e, so we may write it in the view: 

$ E (0) = e4> £ (0). 

Let L £ t (p(u) = E[ip(u £ (t))\u £ (0) = u, as e (0) = x] be the semigroup corresponding 
to the operator L £ . 

Theorem 3.1. The following estimate is true for the remainder (13) of the solution 
of equation (6): 

W(t)\\<e\Wm\^v{eLW\\}, 

where L > 2||(X e )- 1 ||. 

Proof: The solution of equation (14) is: 

f 

**(*)= e[Lf£ e (0) + J L £ t _ s 9 £ {s)ds]. 
o 

For the semigroup we have L £ t = I + L £ J* L^ds, so J Q * L^ds = (L e ) _1 (Lf — I). 
Using Gronwell-Bellman inequality [I], we receive 

t 

||$ £ (t)|| <eLf|||. £ (0)||exp{e | L^(t - s)d S } < eLf ||l» £ (0)|| exp{ e L||0 e ||}, 

o 

where L > 2j|(L £ )- 1 ||. 

Theorem is proved. 
Remark 3.2. For the remainder of asymptotic expansion (5) of the view 

& N+1 (t) :=$ s (t)-& N+1 (t), 

where Q £ N+1 (t) = u^°\t) + ^2^=1 e k (u^ k \t) + w^ k \t)) we have analogical estimate: 

W N+l (t)\\<e N Wm\^v{e N L\\6 £ N \\}, 

where * - L £ <S> £ N+1 := e N 9 £ N . 
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